Voltage rectification by a SQUID ratchet 
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We argue that the phase across an asymmetric dc SQUID threaded by a magnetic flux can 
experience an effective ratchet (periodic and asymmetric) potential. Under an external ac 
current, a rocking ratchet mechanism operates whereby one sign of the time derivative of 
the phase is favored. We show that there exists a range of parameters in which a fixed sign 
(and, in a narrower range, even a fixed value) of the average voltage across the ring occurs, 
regardless of the sign of the external current dc component. 

PACS numbers: 05.40.+j, 74.40. +k, 74.50.+r, 85.25. Cp, 85.25.Dq 
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Although the nonequibrium dynamics of a particle in 
a ratchet potential (i.e., a periodic potential that lacks 
reflection symmetry) has for long been considered a fun- 
damental problem in statistical physics [Q, it has be- 
come the object of more intense attention in recent years, 
because of its newly found relevance in diverse areas of 
physics, chemistry, and biology. A characteristic effect is 
that, when the ratchet is subject to a stationary nonequi- 
librium perturbation, particle motion in one direction is 
favored. Within this context, an important class of dy- 
namical systems is formed by the so-called rocking ratch- 
ets, for which the external perturbation is a time peri- 
odic, uniform force |^^. The effect of dynamically in- 
duced unidirectional motion can overcome the drift effect 
of a small bias that would push the particle into the non 
favored direction. Thus, for not very strong tilts, uphill 
movement is possible provided the ratchet structure is 
conveniently rocked. 

In this letter, we propose a realization of the rock- 
ing ratchet mechanism in a new type of superconducting 
quantum interference device (SQUID) containing a char- 
acteristic asymmetry. The system we propose, depicted 
in Fig. 1, is formed by a ring with two Josephson junc- 
tions in series in one of the arms and only one junction 
in the other arm. We will show that, when the ring is 
threaded by a flux $oxt that is not an integer multiple of 
$0/2 (<I>o = h/2e being the flux quantum), the effective 
potential experienced by the total phase ip across the ring 
displays a ratchet structure. As a consequence, when the 
asymmetric SQUID is 'rocked' by an external ac current 
/(i), one sign of the phase velocity <f is favored. From the 
Josephson voltage-phase relation, we conclude that there 
must be a range of parameters for which a fixed sign of 
the average voltage Vq = h{ip) /2e occurs regardless of 
the sign of the external current dc component Jq. 

We focus on SQUID structures formed by conventional 



Josephson junctions whose phase is a classical variable 
and which can be adequately described by the 'resistively 
shunted junction' model 1^,1). Thus, the phase tp.^ across 
Josephson junction i on the left arm (see Fig. la) obeys 
the equation {i = 1,2): 



Iiit) = Ji sin((/3i) -I- 
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where Ii(t) is the current through the left arm, and Ri, 
Ci, and Ji are the resistance, capacitance, and critical 
current of junction i. For simplicity, we assume here that 
the two junctions in series are identical, and will comment 
later on the case of slightly dissimilar junctions. We take: 
Ci = C2 = 2Cu i?i = i?2 = i?;/2, and Ji = J2 = Ji- The 
total voltage drop across the two junctions is y = V1+V2, 
where Vj, = {h/2e)ipi. 
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FIG. 1. (a) Schematic picture of an asymmetric SQUID 
with three junctions threaded by an external flux, (b) Repre- 
sentation of the equivalent circuit: the two junctions in series 
of the left branch behave like a single junction with ip replaced 
by p/2. 

If ipiit) is a solution for the first junction, then ip2{t) — 
ipi{t) = ipi{t)/2 is also a solution for the second junction 
M. This implies V = (pih/2e, with ipi satisfying the 
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equation: 
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Hence, a series of two identical Josephson junctions can 
be described by the same equation as a single junction, 
with the only difference, that in the sine-function the ar- 
gument Lp/2 occcurs js). This is a most important feature 
to build the ratchet-like structure. On the right arm, the 
phase across the single junction obeys an equation that 
reads as in Eq. (1) with the labels I and i replaced by 
r. In the following we assume that each Josephson link 
operates in the overdamped limit, {2e / K) J ^R^C a -C 1 
(a = /, r), so that the capacitive terms can be neglected 
in Eq. (2) as well as in its right arm counterpart 
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FIG. 2. The ratchet potential (solid) given after Eq. (|) 
of the text, which governs the behaviour of the three junction 
SQUID (of. Fig.|ll), is compared to the sine potential (dashed) 
V(x) = l.lsin(j;). 

The total current through the SQUID is /(<) = //(f) + 
/r(i), and the equivalent circuit 10| is shown in Fig. 
lb. In the limit where the contributions L; and to 
the total loop inductance L — Li + Lr are such that 
\LI\ ^ $0 the total flux $ is approximately the external 
flux $ext- Then, integration of the gauge invariant phase 
around the loop yields tpi — ipr — — (yScxt + 27rn, with 
(/?oxt = 27r$cxt/^o- Thus, we see that the application of 
an external flux provides us with an externally tunable 
relation between ipi and ipr, which, in combination with 
the freedom to choose the ratio Ji/Jr, allows us to select 
the shape of the potential experienced by the phase f = 

We may include the effect of temperature T by adding 
Nyquist noise. We assume Gaussian white noise 'ri{t) of 
zero average and correlation {'r]{t)ri{0)) = {2kBT/ R)S{t), 
so that the phase satisfies the equation 

-^ip = - J; sin (^^^ - Jr sm{ip + ipa^t) + I{t) + ?y(t), 

(3) 

where we have set Ri = R,. = R. The resulting Fokker- 
Planck equation is numerically solved by a matrix con- 
tinued fraction method The total dc voltage across 
the SQUID is given by Uq = {^/2e){ipi) -f L/(/,) = 
{h/2e) {ipr)+Lr{Ir), whcrc (. . .) stands for time and noise 
average. Since currents must remain bounded, one has 
(/;) = (/^) = 0, which leads to Vb = i^/2e) (ip). 
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FIG. 3. The d.c. current-voltage characteristics for the 
SQUID in Fig. § is shown for an adiabatically slow {uj = 0.01) 
a.c. contribution of amplitude ^4 = 1 at different noise 
strengths D: The ratchet potential with noise strengths 
D = (solid) , D = 0.01 (dashed) and D = 0.5 (dotted) 
is compared to the sine potential of Fig. 2 with D — 
(dashed-dotted) and D — 0.01 (dashed-double dotted). Inset: 
Global view of the same I-V curves for the ratchet potential 
with D = 0.01 (dashed) and D = 0.5 (dotted). 

Next we feed the circuit with a current I{t) = Iq + 
Ii sm{Qt). In terms of the dimensionless quantities x = 
(ifi + 7r)/2, T = {eRJi/2h)t, s = Jr/Ji, F = Iq/Ji, A = 
h/Ji, UJ = 2hn/eRJi, and D = eksT/hJi, Eq. (|) reads 

^ = -^U{x)+F + Asm{cjT)+ar), (4) 

where U{x) = — [sin(a;) -I- (s/2)sin(2x -I- fext — t''/2)] is 
the effective potential and ^(r) is Gaussian noise with 
{^{t)£,{Q)) — 2DS{t). The average voltage is now given 
by Vq = {JiR/2) (dx/dr). Setting, for instance, 5 = 1/2 
and ipcxt = 7r/2, U{x) adopts the form of a ratchet po- 
tential with period 27r, as shown in Fig. 2. 

We can expect the ratchet structure arising from the 
combination of asymmetry and ^oxt — <I'o/4 to have 
major consequences on the device properties. In Fig. ^ 
we show the dc current-voltage characteristics for a low 
ac frequency w = 0.01 and A — 1. The resulting dc volt- 
age for the ratchet potential is compared to that obtained 
for a symmetric potential with the same barrier height. 
Clearly, the main effect of the ratchet shape of the poten- 
tial is that of shifting the dc current- voltage characteris- 
tics towards more negative values of Iq or, equivalently, 
towards more positive values of Vq. Within the present 
sign convention, we shall refer to this displacement as a 
shift in the 'right' direction, because in it {(p) has the sign 
that is generally (although not exclusively ^ ; see below) 
favored. 

In the deterministic case, the velocity (dx/dr) is almost 
quantized at values nu, n = 0, ±1, corresponding to 
solutions x{t) that are 'locked' into the phase of the driv- 
ing force [Q. For a symmetric potential, these plateaus 
in the voltage correspond to standard Shapiro 
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FIG. 4. Same as in Fig. ^ for frequency uj — 0.3 and 
A = 1.7, for D = (solid), 0.01 (dashed), and 0.5 (dotted). 
Inset: magnified picture showing steps at fractional values of 
Lj in the D = curve. 



steps [|6yi2|-|14|. In both the symmetric ||15|] and asym- 
metric cases, a small amount of noise {D — 0.01) suf- 
fices to wipe out the structure of steps. However, the 
same weak noise does not destroy the ratchet-induced 
shift in the dc I~V characteristics. As shown in Fig. 3, 
a stronger noise intensity is needed to appreciably re- 
duce the ratchet effect (note that it still persists with D 
as large as 0.5) and to lead the system towards a con- 
ventional behavior in which Vq oc Iq. The same trend 
towards Ohmic response is already shown for weak noise 
if Iq is large enough. 

It is remarkable that the ratchet effect can be displayed 
so clearly at frequencies as low as uj = 0.01. The adia- 
batic limit (uj —>■ 0) can actually be studied analytically. 
For D = 0, one obtains Vq = {JiR/2) /^^ de/t{0), where 
t{e) = dx/[F + Asin(6l) - U'{x)]. For = 0, Vb is 
guaranteed to be zero only if U{x) is symmetric, since 
then U'{x) and thus t{d) must be odd functions. On the 
contrary, if U{x) is not symmetric for any choice of ori- 
gin (ratchet potential), then one generally has Vq 7^ 
with Jq = 0. For a given amplitude, the ratchet behav- 
ior tends to disappear as the frequency grows. On the 
other hand, for a given frequency, there is an optimum 
amplitude that maximizes the ratchet effect 

In Fig. 4, we show Vq as a function of /q for lu = 0.3 
and A = 1.7. In the absence of noise, steps at half-integer 
multiples of w can be clearly observed. In the inset of Fig. 
4, additional steps can be observed at w/m. They are also 
present for uj = 0.01, although they cannot be resolved 
in the scale of Fig. 3. We note that these noninteger 
steps are not due to the ratchet structure itself but to 
the deviation of U (x) from a simple sin(a;) law, which is 
the sole case for which steps lie only at integer values nuj 
p6[ . Under weak noise {D — 0.01), the fractional Shapiro 
steps dissapear, but the structure of integer plateaus is 
still somewhat preserved. More intense noise {D = 0.5) 
destroys the voltage quantization totally and, like in the 



adiabatic case, considerably reduces the ratchet effect. 

So far we have assumed that ac current sources are 
applied to the device. It is interesting to analyse what 
happens when a voltage source of the type V{t) = Vq + 
Vi sm{nt) is applied instead. Then the phase evolves as 



x{t) = a;o + {x)t — {A/uj) cos(cjt). 



(5) 



where (i) /Vq = A/Vi = 2/JiR. Inserting (5) in Eq. (3) 
and averaging over time, one obtains that, for (x) = nw 
or {x) = {2n + l)w/2, a continous interval of dc current 
values F is possible. For ipcxt — one obtains, respec- 
tively, 

F — nuj + Jn 1 — I cos(xi) — sJ2n ( — ) cos(2a;i), 
\ujj \uj J 

2n+l f2A\ ^ ^ 

F= — - — UJ - sJ2n+i \— j sm[2xi), (6) 

where Jn{z) is the n-th order Bessel function The 
finite range of F values spanning a voltage plateau for 
a given value of n is obtained by letting xi take any 
real value. After comparing the structure of plateaus 
predicted by Eq. (6) with that obtained numerically for 
the case of current sources, we have found that, as in the 
symmetric case , similar results are obtained for ui ^ 
1 01 A ^ 1, provided that A/uj < 1. Inspection of Eq. 
(6) shows that the resulting structure of steps, although 
not entirely symmetric, does not exhibit a proper ratchet 
effect in any range of parameters, since there is always a 
Jo = 0, Vb = solution. This can be proved by noting 
that the last two terms in upper Eq. (6) cancel for certain 
values of xi. 

Going back to Figs. 3 and 4, we notice the remark- 
able property that there is a finite range of Iq values in 
which the sign of the average voltage is independent of 
the average external current |l^. For a narrower range 
of parameters (see, e.g., in Fig. 4), and in the absence 
of noise, it is possible to obtain, not only the same sign, 
but also the same value of Vq, regardless of the value and 
sign of Iq. Therefore, we conclude that the asymmetric 
SQUID we propose here can be used as a device for volt- 
age rectification. From the curves presented here, we note 
that this mechanism of voltage rectification will operate 
more efficiently at low frequencies and for not too small 
ac amplitudes jl^ . On the other hand, the analysis given 
in the preceding paragraph indicates that, under the ef- 
fect of an external ac voltage source, the SQUID of Fig. 
1 could not yield current rectification. 

In our analysis, we have assumed for convenience that 
certain ideal relations between the parameters of the dif- 
ferent junctions are satisfied. One may wonder whether 
the physical effects we have discussed may be affected 
by minor deviations from those specific values, especially 
when the two junctions in series are not identical and 
the simple relation (pi{t) = Lp2{t) cannot be always valid. 
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For the case of zero noise, analytical considerations sug- 
gest that a weakened ratchet effect and a structure of 
shorter steps will remain. We have performed a numeri- 
cal check by treating ipi and ip2 as independent variables. 
For LU — 0.01 and A = 1 and 1.7, and assuming dif- 
ferences of order between one and ten percent (namely, 
R2/R1 = J2/J1 = i?r/2i?i = 1.01 and 1.1), we find 
that the dc voltage at zero current bias decreases within 
five to thirty percent and that the voltage plateaus are 
shortened by about one half. These results underline the 
robustness of the predicted physical behavior (in partic- 
ular the ratchet effect) against small deviations from the 
ideal structure. 

For a typical SQUID, the inductance can be L ~ 10^^° 
Henrys jg]. Thus, currents ^ 10~^ A are required for the 
condition $ ~ $oxt to be satisfied. For typical tunnel 
junctions, the overdamped limit cannot be achieved un- 
less one operates very close to the critical temperature 
[^. Operation in a wider range of temperatures could, 
however, be achieved by adding shunts of sufficiently low 
resistance. For J; = 10"^ A, R — 1 Ohm, the 'units' of 
temperature, frequency and voltage are 48 K, 125 MHz, 
and 0.5 fiY. From our numerical results, we may con- 
clude, for instance, that for T = 0.48 K and fl — 37.5 
MHz the dc voltage is Vb ^ O.lfiV at zero dc current. 

In conclusion, we have demonstrated the feasibility of a 
novel effect in the dynamics of the phase across an asym- 
metric SQUID threaded by a magnetic flux. The ratchet 
structure of the effective potential experienced by the 
phase through the ring favors one sign of its time deriva- 
tive. Under an oscillating current source, the dc current- 
voltage characteristics present striking properties such as 
displaced Shapiro steps and the possibility of having a fi- 
nite dc voltage with a zero dc current, and viceversa. 
Within a certain range of parameters, the same sign, and 
even the same value, of the dc voltage can be obtained 
regardless of the sign of the external dc current. This 
mechanism of voltage rectification has been shown to be 
robust in the presence of moderate noise and of small 
deviations of the junction parameters from the proposed 
ideal behavior. Estimates for a single SQUID suggest 
that the predicted ratchet-induced voltage shift is indeed 
measurable. The effect could be conveniently amplified 
by placing many similar devices in series. 
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